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Abstract. We present the first broadband X-ray observations of four Giga-Hertz Peaked Spectrum (GPS) radio 
galaxies at redshift ^ 1 performed by Chandra and XMM-Newton. These observations more than double the 
number of members of this class with measured spectra in hard (E > 2 keV) X-rays. All sources were detected. 
Their radio-to-X-ray spectral energy distributions are similar, except for PKS 0941-080, which is X-ray under- 
luminous by about two orders of magnitude. The comparison between the full sample of GPS galaxies with 
measurements in hard X-rays and a control sample of radio galaxies rules out intrinsic X-ray weakness as causing 
a lower detection rate of GPS sources in X-ray surveys. Four out of seven GPS galaxies exhibit high X-ray column 
densities, whereas for the remaining three this measurement is hampered by the poor spectral statistics. Bearing 
in mind the low number statistics in both the GPS and the control sample, the average column density measured 
in GPS galaxies is larger than in FR I or Broad Line Region FR II radio galaxies, but consistent with that 
measured in High- Excitation FR II galaxies. This leads to a location the absorbing gas in an obscuring "torus" , 
which prevents us from observing the nuclear region along lines-of-sight perpendicular to the radio axis. This 
interpretation is supported by the discovery of rapid (timescale ~ 10 3 s) X-ray variability in the GPS galaxy 
COINSJ0029+3456, and by an almost order-of-magnitude difference between the HI column density measured in 
radio and X-rays in PKS0500+019. 

Key words. Galaxies: jets - Galaxies: active - Quasars:individual:COINSJ0029+3456, PKS0500+019, PKS0941- 
080, PKS2128+048 - X-ray: galaxies 



1. Introduction 

GHz Peaked Spectrum (GPS) radio sources are extremely 
powerful (L ra( uo *S 10 45 erg sec -1 ), very compact (10- 
100 mas, 10-1000 parsecs) sources characterized by a sim- 
ple convex radio spectrum peaking near 1 GHz (O'Dea 
1998 and references therein). GPS sources make up about 
10% of the sources selected at frequencies around 5 GHz. 
They often exhibit symmetric structures on the parsec 
scale, reminiscent of those present in extended radio galax- 
ies on much larger scales. About half of the GPS sources 
are classified as "GPS galaxies". 

GPS sources - together with Compact Steep Spectrum 
(CSS) sources l|0'Dea et al. 1998jl - may play an impor- 
tant role in the history and evolution of the radio uni- 
verse. It was originally suggested that they may repre- 
sent radio galaxies in the early stage of their life (typical 
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ages < 10 4 years; Phillips & Mutel 1982; Carvalho 1985; 
De Young 1993; Fanti et al. 1995; Readhead et al 1996; 
Murgia 2003). This possibility was recently supported by 
the detection of micro arc-seconds hotspot proper motions 
in a sample of ~10 GPS galaxies (Poladitis & Conway 
2003). However, evidence for weak ~10 kpc scale radio 
emission in 5-10% of all known GPS (Stanghellini et al. 
1998 and references therein) suggested that the radio ac- 
tivity may be either recurrent (Baum et al. 1990), or 
smothered by the inflow of "fresh" matter after merging 
or interaction with a gas rich companion. Alternatively, as 
originally suggested by Gopal-Krishma & Wiita (1991), 
GPS sources could remain compact during their radia- 
tive lifetime because interaction with dense circumnuclear 
matter impedes their full growth. Semi-analytical models 
(Carvalho 1994; 1998) and hydrodynamic simulations (De 
Young 1993) have indeed shown that InterStellar Medium 
(ISM) average densities in the range 1-10 cm -3 could be 
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very efficient in preventing the development of large scale 
radio structures. 

Measurements in X-rays could provide fundamental 
clues on the nature of GPS sources and consequently on 
the origin and evolution of the radio power in the universe. 
A hot (T ~ 10 6 ~ 7 K) and tenuous ISM phase with n ~ 
1 cm -3 should be a copious source of soft X-rays. However, 
the recent Chandra discovery of the X-ray cluster associ- 
ated with the CSS quasar 3C 186 rules out the confine- 
ment of the radio source by external medium in this quasar 
( Sicmigi nowsha et al. 2005al ). Alternatively, a distribution 
of cold clouds would imprint an X-ray photoelectric cutoff 
on the soft X-ray spectrum. It has been noted that GPS 
sources are rather elusive in X-rays. In the ROSAT All 
Sky Survey GPS/CSS quasars exhibited a detection rate 
three times lower than radio-loud quasars of comparable 
power ( Bake r~et al. 1995|l . Whether this low detection rate 
is due to intrinsic weakness or to obscuration of the ac- 
tive nucleus is still largely unknown. Prior to the launch of 
Chandra and XMM-Newton, high-quality CCD-resolution 
hard (E > 2 keV) X-ray spectra were available for two 
GPS galaxies o nly: NGC 1052 llOuainazzi et al. 20n71f) 1 
and 4C+12.50 <|Q'Dea et al. 2000|1 . Both exhibit a high 
X-ray obscuring column density (Nh > 10 22 cm -2 ). 

In order to gain insight into the X-ray under- 
luminosity of GPS galaxies , we have performed Chandra 
and XMM-Newton observations of 5 GPS galaxies at 
z £ 1. This paper presents the results of this program. 
The observations discussed here represent the first hard 
X-ray measurements of the target sources. They more 
than double the number of GPS galaxies with hard X-ray 
measurements. The results of the observation of Mkn 668 
(OQ+208), performed in the framework of our program, 
were presented by Guainazzi et al. (2004). This observa- 
tion allowed the discovery of the first radio-loud Compton- 
thick {Nh > ai 1 ) AGN. 

Unless otherwise specified: energies are quoted in the 
source reference frame; uncertainties on the spectral pa- 
rameters are at the 90% confidence level for one interest- 
ing parameter; upper limits are also at the 90% confidence 
level; other uncertainties are at the la level; in the calcu- 
lation of the luminosities, we adopted a Hubble constant 
of 70 km s _1 Mpc -1 l|Bennett et al. 20 03): data points 
in the spectral plots correspond to a signal-to-noise > 3. 
We have used throughout the paper the Gehrels algorithm 
IjGehrels 1986|l to estimate the statistical uncertainties as- 
sociated with Poissonian-distributed measurements. 

2. The sample 

In this section we present a short description of the 
radio properties of the sources discussed in this paper, 
supporting their classification as "GPS" galaxies. We 
require both a convex radio spectrum and a symmetric 
radio morphology to classify a radio galaxy as GPS. 

1 see Sect. 5.2 for a discussion on the GPS classification of 
this radio galaxy 



COINS J0029+3456 (B 0026+3 46) is a galaxy 
with m r =21.0 2 at z=0.517 l|Snellen et al. 19961 
IZensus et al. 2002f . It has a broad radio spectrum, which 
is rather flat at high frequencies l|Kiihr et al. 1981|) . It 
exhibits a compact symmetric morphology with an overall 
size of 35 mas (190 pc). 

PKS 0500+019 (J0 503+0203) is a ra dio galaxy with 
m.^21.0 at z=0.583 l|Carilli et al. 1 998 ) with a Compact 
Symmetric Object (CSO) morphology which extends for 
about 15 mas (84 pc). Stickel et al. (1996) reported an 
emission line at a higher redshift, therefore they consider 
the GPS radio source associated with a background 
quasar and the galaxy as an intervening object. De Vries 
et al. (2000), however, did not find any trace of this 
emission line in their optical spectra. 

PKS 0941-080 is a radio galaxy at z = 0.228 with a 
clear convex radio spectrum and a rather clear CSO 
morphology (Dallacasaet al. 1997). 

PKS 2128 +048 (J2130+0502) is a galaxy with m r =23.3 
at z=0.99 l|Snellen et al. 1996|> with a clear GPS radio 
spectrum and a rather clear CSO morphology of 35 mas 
(220 pc) UStanghellini et al. 1997| 2001). 

Table ^ reports the list of observations discussed in 
this paper. The March 2004 observation of PKS0500+019 
was heavily affected by high-particle background, which 
reduced the usable exposure time to ~2 ks. The target 
was therefore successfully re-observed in August. 

3. XMM-Newton results 

The XMM-Newton data were reduced with SAS v6.0 
l|Gabriel et al. 2 0031. using the most updated calibration 
files (January 2005). In this paper, only data from the 
EPIC cameras (MOS; Turner et al. 2001; pn, Striider et 
al. 2001) will be discussed. Event lists from the two MOS 
cameras were merged before accumulation of any scien- 
tific products. Single to double (quadruple) events were 
used to accumulate pn (MOS) spectra. High-background 
particle flares were removed by applying fixed thresh- 
olds on the single-event, E > 10 keV, At = 10 s light 
curves. These thresholds, as well as the radius of the 
source scientific product circular extraction regions, were 
optimized to maximize the signal-to-noise ratio, and are 
listed in Tabled Background scientific products were ex- 
tracted from annuli around the source for the MOS, and 
from circular regions in the same chip for the pn, at the 
same height in detector coordinate as the source location. 
Spectra were binned in order to oversample the intrinsic 
instrumental energy resolution by a factor >3, and to have 
at least 25 counts in each background-subtracted spectral 
channel. This ensures that the % 2 statistics can be used 
to evaluate the quality of the spectral fitting, pn (MOS) 
spectra were fitted in the 0.35-15 keV (0.5-10 keV) spec- 
tral range. 

2 m r and rrii are the magnitudes in the r and i band, respec- 
tively, according to the Gunn system. 
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Table 1. Log of the observations discussed in this paper. "Threshold" indicates the maximum background count rate 
used in the generation of appropriate time intervals for the extraction of spectra. The "Radius" refers to the source 
extraction region 



Source 


z 


Nh,Go.i 


Start date 


-L exp 


Threshold (s _1 ) 


Radius (") a 






(10 20 cm" 2 ) 




(ks) 






COINS J0029+3456 


0.571 


5.6 


08-Jan-2004 


12.2/10.8 


0.25/1.0 


35/40 


PKS 0500+019 (1) 


0.584 


8.3 


09-Mar-2004 


2.0/1.9 


2.0/1.0 


20/30 


PKS 0500+019 (2) 






18-Aug-2004 


9.6/7.2 


0.0/0.0 


20/30 


PKS 0941-080 


0.228 


3.7 


26-Mar-2002 


5.3 




1.4 


PKS 2128+048 


0.990 


5.2 


ll-Oct-2002 


5.7 




2.4 



"for XMM-Newton MOS/pn and Chandra ACTS, respectively 



Table 2. GPS sample sources X-ray coordinates, "d" 
indicates the distance between the X-ray and the ra- 
dio coordinates, the latter extracted from the NED cat- 
alogue. "CTS" are net counts in the 0.5-10 keV and 
0.1-11 keV energy bands for the XMM-Newton EPIC 
and the Chandra ACIS, respectively. The values for 
PKS 0500+019 refer to the second observation only 



COINSJ0029 +3456 Bin time: 1024. s 



Source 


RA 


Dec 


d(") 


CTS 




(J2000) 


(J2000) 






0029 


00 h 29 m 14M 


+34°56'32" 


1.7 


1030 ± 40 


0500 


05 h 03 m 20 s .9 


+02°03'05" 


4.5 


1350 ± 40 


0941 


09 h 43 m 36 s .9 a 


-08°19'31" a 




10 + 3 


2128 


21 h 30 m 32 s .8 


+05°02'08" 


1.1 


92 + 10 



I fixed to the radio position 



3.1. Source identification 

In Table El we list the coordinates of the nominal targets 
in our XMM-Newton observations. They were determined 
on the 0.5-10 keV image extracted from merged event 
lists of the three EPIC cameras. Typical statistical er- 
rors are ~1", whereas typical systematic uncertainties 
on the absolute attitude reconstruction are ~1.5". For 
PKS0500+019 (2 nd observation) the X-ray position is 
shifted by about 4.5" in RA with respect to the position 
of the radio source IjCondon et al. 19981 IMa et al. 1998|l . 
corresponding to about 25 kpc at the source redshift. 
This difference exceeds the expected measurement un- 
certainties. The XMM-Newton position agrees well with 
that of the closest source in the ROSAT PSPC catalog 
jVoges et al. 2000| | . No optical counterpart exists for ei- 
ther the X-ray or the the radio position. The paucity of 
X-ray/optical coincidences in the XMM-Newton field pre- 
vented us from performing an independent calibration of 
the absolute astrometry in this observation. Nonetheless, a 
comparable shift in RA is observed between several XMM- 
Newton sources and the closest Digital Sky Survey plate 
source. Although this evidence does not constitute proof, 
we will consider in the following such a shift as due to 
residual inaccuracies in the attitude reconstruction. 



Net EPIC 



t 3 



~T~ 



t 1 — i r~ 



i 1 — I 1 — i 1 — r~ 



i 1 1 1 r 



j I i i i i I i_ 



_i I i i i i I 



o 



5000 



10 4 



1.5xl0 4 



2xl0 4 



Time (s) 

Start Time 13012 19:03:40:037 Stop Time 13013 0:27:56:037 

Fig. 1. 0.5-10 keV EPIC (pn and MOS) light curve dur- 
ing the XMM-Newton observation of COINS J0029+3456. 
The binning time is At = 1024 s. Times are in seconds af- 
ter the observation start. 

3.2. Timing analysis in COINS J0029+3456 

The light curves observed during the XMM-Newton ob- 
servations are consistent with being constant in all energy 
ranges, except for COINS J0029+3456, where variability 
of a factor of ~3 within a few thousand seconds was ob- 
served (Fig.^). By fitting consecutive segments of the light 
curve with a linear function, we estimate a typical lumi- 
nosity change rate of (5.0 ± 2.5) x 10 43 erg s -1 hour -1 . 
We find only marginal evidence that the variability is pri- 
marily due to hard X-rays. The reduced x 2 resulting from 
fitting with a constant line the At = 1024 s binned light 
curves is 0.86 and 1.27 in the 0.2-1 and 1-10 keV energy 
band, respectively. Consequently, no significant variability 
of the hardness ratio is associated with these flux changes. 
We will therefore discuss the time-averaged spectrum of 
this source in Sect. 3.3. COINS J0029+3456 also shows a 
remarkable historical variability in the soft X-ray band. 
The 1 keV flux density decreased from 100 ± 30 nJy 
(Einste in/IPC, 1980) and 63 ± 15 nJy (ROSAT/PSPC, 
1992) (Koll gaard et al. 1995| our re-analysis of archival 
data) to 29±7 nJy as measured about 10 years later 
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(January 2004) by XMM-Newton/EPIC. It is impossible 
to tell whether a change of the column density of the in- 
tervening absorber contributes to this variability, due to 
the lack of statistical quality and energy resolution of the 
measurements prior to XMM-Newton. 

3.3. Spectral analysis 

Spectra of the EPIC cameras have been simultaneously 
fitted, allowing only a free cross-normalization factor be- 
tween them. We have applied a "baseline model", con- 
stituted by a power- law modified by photoelectric absorp- 
tion. We have modeled the absorption simultaneously with 
two components: a) the intervening gas along the line of 
sight in our Galaxy, with a column density Nn,Gai held 
fixed to the value determined by the 21 cm radio maps 
( |Dickey fc Lockman 1990] ) and; b) a "local" absorber lo- 
cated at the same redshift as the systemic velocity of each 
galaxy; its column density, Nh, was left as a free param- 
eter in the fit. All the continuum parameters were left 
free in the fit as well. This baseline model yields a good 
fit for all observations. The spectral parameters derived 
from the two observations of PKS 0500+019 are formally 
consistent within the statistical uncertainties, except for a 
~40% larger average flux in the later one. We will consider 
in the following the results obtained during the second ob- 
servation of PKS 0500+019, as their associated statistical 
uncertainties are smaller. 

Table summarizes the spectral parameters of the 
XMM-Newton GPS sample presented in this paper. The 
spectra and residuals against the best fit models are shown 
in Fig. El 



4. Chandra results 

Two GPS galaxies, PKS 0941-080 and PKS 2128+048, 
were observed with the Chandra Advanced CCD Imaging 
Spectrometer (ACIS-S, [Weisskopf et al. 2002) ) as part 
of a program to observe a sample of GPS sources 
(Siemiginowska et al., 2005b). The short 5 ksec observa- 
tions were mainly scheduled to measure the X-ray flux. 
The sources were located ~ 35" from the default aim-point 
position (to avoid node boundaries) on the ACIS-S back- 
side illuminated chip S3 (Proposer's Observatory Guide 
3 ). The 1/8 sub-array CCD readout mode of one CCD 
only was used resulting in a 0.441 sec frame readout time. 

The X-ray data analysis was performed in CIAO 
3.1 4 with the calibration files from the CALDB 2.28 
data base. Note that the ACIS-S contamination file 
acisD1999-08-13contamN0003 . f its which incorporates 
effects of the ACIS-S contamination into the effective area 
was included in our analysis. The spectra were rebinned 
according to the same criteria employed for XMM-Newton 
(sec Sect. 3). 



2 10 - 




0.0 



0.5 1.0 1.5 

Radio spectral index (>'>v breal< ) 



2.0 



3 http:/ /asc. harvard.edu/proposer/POG/index. html 

4 http://cxc.harvard.edu/ciao/ 



Fig. 4. 5 keV to 5 GHz flux density ratio as a function of 
the radio spectral index above the break frequency. The 
slope of a linear least squares fit is (—1.2 ± 0.6) x 10~ 8 . 
The slope does not significantly change if the data point 
corresponding to PKS 0941+080 is excluded from the fit. 



The two galaxies were detected by Chandra with the 
count rates listed in Table For PKS 0941-080 the to- 
tal number of net counts (~10) is too low to perform a 
true spectral fitting. In PKS 2128+048 a simple power- 
law model with photoelectric absorption yields a good fit. 
A column density of Nh — 3 x 10 21 cm -2 is measured. 
However this measurement is still consistent with if the 
large statistical uncertainties are taken into account. 



5. Discussion 

5.1. Spectral Energy Distributions 

In Fig. |31 we show the radio-to-X-ray Spectral Energy 
Distributions (SEDs) for the GPS galaxies discussed in 
this paper, as well as for 4C+ 12.50 IjO'Dea et al. 2000jl 
and Mkn 668 IjGuainazzi et al. 20040 . The 5 GHz to 5 keV 
flux density ratios of 5 out of 6 objects are within a 
factor of 15. On the other hand, in PKS 0941+080 the 
5 GHz to 5 keV flux density ratio is about two orders- 
of-magnitude lower than the average for the rest of the 
sample. Whether this is due to intrinsic weakness or ex- 
treme (i.e. Compton-thick) obscuration is impossible to 
tell with the current data, and would require very deep 
X-ray exposures to be finally elucidated. There is no cor- 
relation between the X-ray to radio flux density ratio and 
the radio break frequency; however, an anti-correlation at 
the 2cr level is found between the X-ray to radio flux den- 
sity ratio and the radio spectral index above the frequency 
break ( |Stanghellini et al. 19981 see Fig. 00. 

Mkn 668 is a peculiar object in our small sample. 
Not only does it exhibit the largest measured obscura- 
tion among the GPS galaxies for which hard X-ray mea- 
surements are available (being one of the very few known 
radio- loud Compton-thick AGN), it also shows a promi- 
nent soft X-ray excess above the extrapolation of the ob- 
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Table 3. Spectral results on the XMM-Newton GPS sample described in this paper. The model employed to fit the 
data is a photoelectrically absorbed power-law. The obscuring gas with column density Nh is assumed to be located 
at the systemic velocity of the galaxy. T is the photon index of the power-law, F the observed flux, L the intrinsic 
luminosity, corrected for absorption. 



Source Nh f F 1 I? x 2 /v 

(10 21 cm- 2 ) 

COINSJ 0029+3456 10±i 1.43±!J;S 2/7 2.3 ±0.2 38.9/34 

PKS 0500+019 (1) 10+p 1 2.0±8:I 3.9 3.8 + 0.9 10.3/9 

PKS 0500+019 (2) 5±| 1.62+g:^ 5.6 5.0 + 0.6 33.1/44 

PKS 0941-080 ... 2 C 0.07 0.009 + 0.008 ... d 

PKS 2128+048 3.0±§;o 1.5±o:? 1.6 4.4 + 1.1 1.7/4 



a in units of 10 13 erg cm 2 s 1 in the 0.5-10 keV energy band (observer's frame) 

6 in units of 10 44 erg s _1 in the 2-10 keV band (source frame), corrected for absorption 

c fixed 

d not enough independent channels for a spectral fit to be possible 



COINSJ0029+3456 PKS0500+019 (2) 




0.5 1 2 5 10 0.5 1 2 5 10 

Energy (keV) Energy (keV) 



Fig. 2. Spectra {upper panels) and residuals in units of standard deviations [lower panels) when the best fit model of 
Tablets applied to the XMM-Newton/EPIC spectra of COINS J0029+3456 and PKS 0500+019. 



scurcd primary AGN continuum. None of the other GPS 
exhibits this component in the XMM-Newton or Chandra 
spectra. Its origin is still to be fully elucidated, the most 
likely possibilities being scattering of the AGN radiation 
or inverse- Compton recoil of the large far infrared emis- 
sion (fMazzarel la et al. 19911 IGuainazzi et al. 2 004 ) . 

5.2. Comparison with a control sample of "normal" 
radio galaxies 

In this Section we compare the spectral properties 
of the full sample of GPS galaxies for which hard 
X-ray measurements are available with a control sam- 
ple of Radio Galaxies (RGs). The former includes, 
with the objects discussed in this paper: 4C+12.50 
IjO'Dea et al. 20000, 4C+55.16 Ijlwasawa et al. 1999|) . 
OQ+208 l|Guainazzi et al. 20040 . We do not in- 
clude the radio galaxy NGC 1052, despite it be- 
ing sometimes classified as a "recurrent" GPS 
source QTingay et al. 2003| |Edwards fc Tringay 20040 . 
NGC 1052 has a 2-sided jet morphology at the 
mas scale, without evidence of micro hot spots 



IjKadler et al. 20041 |Satypal et al. 20 040 . The hot spots 
are instead present at large scales, implying that the the 
jets are a continuous structure, which ends at the kpc 
scale. 

The control sample is taken from a compilation of 
radio-loud AGN X-ray spectra discussed by Sambruna et 
al. (ASCA, 1999; S99 hereafter), Donato et al. (Chandra, 
2004; D04), and Fiocchi et al. (BeppoSAX, 2005; F05). 
S99 provides in the redshift range of the GPS sample a 
not statistically complete sample of 15 radio galaxies of 
mixed FR I and FR II types, and optical (broad/narrow 
line) properties. To this, F05 add five more FR II and D04 
three more FR I galaxies, again of mixed optical types. 
This comparison aims to address the following basic ques- 
tions: a) are GPS galaxies intrinsically X-ray weaker than 
"normal" radio galaxies?; b) are GPS galaxies more ob- 
scured than "normal" radio galaxies? 

5.2.1. Are GPS galaxies intrinsically X-ray weaker than 
"normal" radio galaxies? 

No. In Fig. we show the location of the GPS sample and 
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Fig. 3. Radio-to- X-ray Spectral Energy Distributions for the objects of our sample, 4C+12.50 (JO'Dea et al. 2000) 
and Mkn 668 ( Guainazzi et al. 2004). Radio measurements represent not simultaneous observations compiled from the 
NASA Extragalactic Database (NED). 



of the control radio galaxy sample in the 5 GHz versus 
2-10 keV luminosity plane (we stress that the X-ray lu- 
minosity is intrinsic, i.e. corrected for absorption). GPS 
galaxies occupy the same region of the Lx versus L ra( n 
plane as "normal" radio galaxies across the whole range 
of X-ray luminosities covered by the sample. The average 
luminosity ratios Rx r are: \og[Rxr{GPS)\ = 0.7 ± 0.3, 
\og[R Xr (RG)} = 0.69 ± 0.17. 



The origin of the X-ray continuum in GPS galaxies is 
still largely unknown. Observed spectral indices are in the 
range 1.4-2.0, therefore consistent with both accretion- 
IjPerola et al. 20021 IRisaliti 2002|) or beamed-dominated 
IjFossati et al. 19 98) systems. The upper limits on the 
equivalent width of the K Q fluorescent iron line (^150 eV) 
are not stringent enough to rule out the former explana- 
tion. Moreover, for none of the objects presented in this 
paper are measurements at 7000A available, which may 
contribute to elucidate the nature of the central core con- 
tinuum through the radio-to-UV correlation in Chiaberge 
et al. (1999). 



5.2.2. Are GPS galaxies more obscured than "normal" 
radio galaxies? 

This may be true. In a recent Chandra survey of FR I RGs 
(D04), only one out of 15 objects exhibit a column density 
~10 22 cm -2 , whereas for 13 of them Nh < 10 213 cm~ 2 
(D04). In Fig. we show the measured column densities 
for the GPS sample and the FR II sample of S99 and 
F05, together with the Nh distribution function for the 
two classes. The apparent anti-correlation between Nh 
and z is most likely due to a selection effect, as farther 
objects with greater obscuration are more easily missed. 
All the measurements of Njj in GPS galaxies lie on the 
upper envelope of the plot. With the caveat of the small 
numbers, the two distributions are different at the 98.3% 
level, according to the Kolmogorov-Smirnov (K-S) test. 
The difference is primarily due to the relative lack of un- 
absorbed GPS galaxies. Only one out of five GPS galaxies 
(< 46%) has Nh £ 10 22 cm" 2 , compared to 75 ± 26% 
RGs. If the X-ray weakness of PKS 0941-080 is due to 
high obscuration, or we include in the sample NGC 1052 
(Nh — 2 x 10 23 cm" 2 . [Guainazzi et al. 20 00) the evidence 
for greater absorption in GPS galaxy is even more signifi- 
cant (see as well Guainazzi et al. 2004). These results are 
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Fig. 6. Left panel: Nh versus z distributions for GPS galaxies (filled circles) and the control sample of FR II radio 
galaxies (empty squares). Right panel: the Njj distribution functions for FR II radio galaxies (top) and GPS galaxies. 
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Fig. 5. 5 GHz versus 2-10 keV for the control sample ra- 
dio galaxies (empty squares) and the GPS galaxies sam- 
ple (filled circles). Empty circles surround the FR I radio 
galaxies data points. 



still sensitive to Poissonian noise. Their confirmation on a 
larger sample is a task we are actively pursuing. 



Neither the control sample nor the GPS sample are, 
by construction, unbiased or complete. However, they are 
reasonably well matched in both radio luminosity and red- 
shift. The K-S probability that the radio luminosity distri- 
butions for GPS and the control sample are drawn from 
the same parent population is 66.5%. The same proba- 
bility for the redshift distribution is somewhat smaller 
(37.4%; cf. Fig. 0). This is potentially more worrying, as 
the GPS redshift distribution median is larger than in the 
control sample, and obscured objects are in principle more 
easily detected at higher redshifts. In order to estimate 
the potential systematic error due to comparing column 
density distributions corresponding to samples with differ- 
ent redshift distributions, we have simulated a large num- 
ber of spectra, using the best-fit models applicable to the 
GPS objects of our sample, but assuming that they are 
distributed in redshift in the same way as the objects of 
the control sample. The difference between the observed 
and the simulated GPS Nh distributions were parameter- 
ized through the statistical uncertainty-weighted average 
of the shift in the column densities, ANh- The shifts are 
-0.19 ± 0.09 dex and -0.08 ± 0.09 dex only, when the 
redshift distributions for FR II and FR I are used. 

Previous studies in the optical and UV revealed a high 
detection rate for FR I RGs ( |Ghiaberge et al. 1999| 2002). 
Their X-ray spectra are generally unobscured (D04) The 
emission of their central compact core has been inter- 
preted as non-thermal synchrotron radiation from the base 
of the jet, suggesting that the we are observing the in- 
nermost regions close to the supermassive black hole. In 
the framework of the grand-unification scenario for FR 
sources flUrry &: Padovani 1995| ) the larger obscuration of 
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Fig. 7. Redshift distribution for the GPS sample (lower 
panel) and the "normal" radio galaxy control sample (up- 
per panel) 

GPS galaxies could be ascribed to their orientation with 
respect to an azimuthally-symmetric absorbing structure 
("torus" hereafter), which prevents us from directly ob- 
serving the central regions in objects observed along lines- 
of-sight perpendicular to the jet axis. 

More puzzling is the evidence that GPS galaxies are 
"more obscured" than FR II objects as well. According 
to a classification originally due to Jackson & Rawlings 
(1997), FR lis can be distinguished among Broad Line 
Objects (BLO), and High- and Low-Excitation galaxies 
(HEG, LEG; |Laing et al. 1994|). BLQs are interpreted as 
unobscured nuclei ij'Varano'^'t al. 2004jl . Optical nuclei of 
LEGs are detected in a comparable fraction as FR I, and 
are therefore more similar to that class (which requires 
some revision of the standard zero-th order unification 
scenario). Finally, 50% of the HEG do not show any op- 
tical counterpart, and should represent the "true" torus- 
obscured radio-loud AGN population. The covering frac- 
tion of this "torus" seems to be redshift-dependent, its 
opening angle widening up as the accretion disk luminos- 
ity increases l|Varano et al. 20041 see as well Page et al. 
2003; Jimenez-Bailon et al. 2005). In our control sample, 
both BLO (5 objects) and LEG (2) FR II RGs exhibit typ- 
ically N H <, 3 x 10 21 cur 2 . The HEG (4) show a much 
larger spread in column densities (Nh & 10 23 cm' 2 ), with 
a {log{N H ^ HEG )) — 2.7 ± 0.6, consistent with that of the 
GPS galaxy sample [(log(N H ,HEG)) > 2.2 ± 0.7]. At face 
value, this implies that the line-of-sight to HEG FR II and 
GPS galaxies AGN intercepts matter of comparable den- 
sity, suggesting a common orientation with respect to the 
obscuring "torus" postulated for the former. No measure- 



ment of the orientation of the radio structure is available 
for the objects of the GPS sample. However, their sym- 
metric morphology suggests that the inclination angle is 
not small. If X-rays originate close to the supermassive 
black hole, and a torus surrounds the central region, our 
line-of-sight will very likely intercept it. 

Most of the GPS host galaxies in our sample exhibit 
a distorted morphology (4C+12.50, PKS 0500+019, 
Mkn 668) and/or double nucleus or close companions 
(all except COIN S J 0029+3456, IPeacock et al. 1981 
Bire tta et al. 19851 |Fugmann fc Meisenheimer"l 988 
Stanghellini et al. 1993| Ide Vries et al. 1995 

Stickel et al. 199611 . High-resolution 2.2/x NICMOS/HST 
images of 4C+12.50 and PKS 0941-080 confirm the 
presence of double nuclei and, in the former case, of 
a highly disturbed morphology (Stanghellini et al., 
submitted; Guainazzi et al., in preparation; see as 
well de Vries et al. 1998, 2000 for ground-based near- 
infrared observations). These reports indicate recent 
mergers. The role of gravitational torques induced 
by galaxy mergers in driving gas inflows toward the 
nucleus has been highlighted by numerical simulations 
dBarnes fc Hernquist 1996) |Taniguchi fc Wada 19961 ) . 
Indeed, interacting AGN pairs seem to invariably host 
X-ray obscured nuclei (see Guainazzi et al. 2005, and 
references therein). Galaxy mergers could therefore 
represent the link between the large amount of gas in 
the GPS nuclear environment and the disturbed and 
multiple optical morphologies of GPS host galaxies 
IJO'Pea etal. 199611. 

Does the gas probed by X-ray measurements have 
enough density to influence the evolution of the parsec- 
scale jet? Guainazzi et al. (2004) demonstrated in the 
cases of Mkn 668 that if the parsec-scale jet had to pierce 
its way through the Compton-thick absorber covering the 
innermost nuclear region, this could lead to an underesti- 
mate by up to two orders of magnitude of the jet dynam- 
ical age as derived from the hotspot recession velocity. 
Even a Compton-thick absorber would not be enough to 
ultimately "choke" the expanding jet. However, measure- 
ments of HI absorbing gas ( Pihlstro m et al. 2003(1 suggest 
a particle density of n e ~30 cm -3 at the external surface 
of the "drilling jet" (consistent with a residual column 
density covering the soft X-ray excess) , which is still con- 
sistent with the ISM density needed to "frustrate" its evo- 
lution. For the GPS galaxies presented in this paper, radio 
HI measurement are available only for PKS0500+019. The 
r ~84 pc separated hot spots are covered by a HI column 
density of ~ 6 x 10 20 cm' 2 , therefore about one order 
of magnitude lower than observed in X-rays. This implies 
that the X-ray absorber is located inside the hot-spots. If 
the jet has to drill its way through it before reaching the 
current separation, the expansion time is: 

t e <, 7.3 x 10 4 L^ 5 0?o 5 years 

HScheuer 19741 ICarvalho 1985fl where L44 is the (un- 
known) luminosity injected in the jet in units of 10 44 and 
fiio is the jet opening angle in units of 10°. The "braking 
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Fig. 8. Nh versus X-ray luminosity for GPS galaxies 
(stars) and GPS/CSS quasars (filled circles). The small 
triangles represent the radio-loud quasars observed by 
ASCA after Reeves & Turner (2000) 

effect" of the X-ray absorber is therefore negligible. Still, 
any residual particle density facing the expanding jet head 
is constrained only to be n e <5 Nh /r = 20 cm~ 3 . We know 
very little of the distribution of the ISM in the innermost 
kpc of this (and other) AGN and the ultimate fate of the 
jet in PKS0500+019 is not known. 

With typical intrinsic X-ray luminosities 
~10 44 erg s _1 , most GPS galaxies of our sample 
have a borderline luminosity between the galaxy and the 
quasar regimes. The X-ray spectral properties and lumi- 
nosities suggest that they may belong to the long-sought 
class of obscured high-luminosity sources predicted by 
theoretical models of the Cosmic X-ray Background 
(CXB; IGilli et al. 200l)l . The cosmological evolut ion of 
GPS sources - which represents a sizable fraction of 
radio-selected extragalactic sources HQ'Dea et al. 1998JI - 
may have an impact on the CXB. In Fig. [5] we compare 
the Nh versus X-ray luminosity in the GPS galaxy sample 
with that measured in a sample of GPS/CSS quasars 
(S99; Siemi gmowska et al. 2003| F05), and in the large 
compilation of spectra of radio-load quasars observed by 
ASCA after Reeves & Turner (2000). In the latter, a 
correlation exists between the obscuring column density 
and the absorption-corrected 2-10 keV luminosity. Given 
the large fraction of upper limits on Nh in the sample, 
we have applied a linear fit log(A#) = a + b x \og(Lx) 
by extending the regression method for left censored 
data described by Schmitt (1985) and Isobe et al. (1986). 
We performed a large number of least square fits on 
a set of Monte-Carlo simulated data derived from the 



experimental points according to the following rules: a) 
each detection was substituted by a random Gaussian 
distribution, whose mean is the best-fit measurement and 
whose standard deviation is its statistical uncertainty; 
b) each upper limit U was substituted by a random 
uniform distribution in the close interval [0,U]. The 
slope obtained with this method - albeit rather shallow 
- is larger than at a confidence level greater than 
ha: b QSO = (3.9 ± 0.7) x 10" 3 . The same fit applied 
on the GPS/CSS sample yields no significant evidence 
for a correlation, even if we include the GPS galaxies 
discussed in this paper: b GPS / css = (— 30 ± 100) x 10~ 3 , 
although the large statistical uncertainties prevent us 
from drawing any firm conclusions. If this different 
behavior is confirmed on a larger sample, this may 
indicate that GPS/CSS quasars are a "mixed bag" of 
obscured/unbeamed and unobscured/beamed objects. 
Beamed structures in high-luminosity quasars may create 
an unobscured line-of-sight by piercing their way through 
dense nuclear environment gas. In statistical terms, this 
would happen more rarely in lower-luminosity galaxies. 
This scenario cannot be investigated with the currently 
available sample, and would require good quality X- 
ray observations on a significantly larger, statistically 
unbiased sample of GPS sources. 

6. Conclusions 

This paper investigates the apparent X-ray under- 
luminosity of GPS sources shown by ROSAT 
( Ba ker et al. 1995|l . by obtaining for the first time 
good quality X-ray spectra on a sizable sample thereof. 
The main conclusions of this study are: 

— The GPS radio-to- X-ray Spectral Energy distributions 
are generally similar. The ratio between the 5 GHz and 
the 5 keV flux density is constrained within a factor 
~15. There is no unique reason for this residual scat- 
ter, which is likely to be due to a combination of X- 
ray obscuration and intrinsic luminosity effects. On the 
other hand, PKS 0941-080 is two orders of magnitude 
X-ray under- luminous with respect to the average of 
the other objects discussed in this paper. Whether this 
is due to high obscuration or intrinsic X-ray weakness 
is impossible to tell with the current data. We there- 
fore refrain from inferring any implications from this 
single outlier 

— The detection of rapid variability during the XMM- 
Newton observation of COINS J0029+3456 puts a 
constraint on the high-energy photon emitting region 
^10 /xpe, supporting an origin at the base of the ex- 
panding jet (where radio synchrotron emission is likely 
to be self-absorbed) or in the accretion disk surround- 
ing the active nucleus. 

— The obscuring screen in PKS0500+019 is located well 
within the radio hot spots, as suggested by the fact 
that the column density determined in X-rays is almost 
one order of magnitude higher than that estimated 
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from radio measurements. The hot spots are therefore 
completely X-ray silent in PKS0500+019. The lack of 
soft excess emission above the obscured X-ray emission 
in COINS J0029+3456 can be interpreted likewise 

— GPS galaxies exhibit a radio-to- X-ray luminosity ra- 
tio comparable to a control sample of "normal" radio 
galaxies. This indicates that it is unlikely that GPS 
galaxies are intrinsically X-ray weaker than their less 
compact counterparts 

— The X-ray obscuring gas column density distribu- 
tion function for GPS galaxies is inconsistent at the 
<;98.7% level with the distribution function of a con- 
trol sample of RGs. The main reason for this differ- 
ence is the comparative paucity of X-ray unobscured 
(N H < 10 21 - 5 cm" 2 ) GPS galaxies 

— The average X-ray column density in our sample is con- 
sistent with that observed in the 4 HEG FR II galaxies 
for which measurements of this observable are avail- 
able. This suggests that X-ray absorption in a GPS 
galaxy is mainly due to orientation with respect to an 
obscuring "torus" 

— There is no compelling evidence for hot (kT ~ 10 6 - 
10 7 K) gas from the X-ray spectra 

The above evidence support a scenario whereby GPS 
radio galaxies are "young" counterparts of large-scale 
FR II radio galaxies. Nh ~ 10 22 cm -2 column densities in 
the innermost kpc around the active nucleus are unlikely 
to substantially affect the evolution of sub-kpc jet struc- 
tures. However, tenuous ISM matter could have a much 
more important role in braking or smothering the full ex- 
pansion of the radio structure. Signatures of this hot gas 
could be revealed by deep high-resolution X-ray Chandra 
observations. 

Although the GPS sample presented in this paper 
represents the largest compilation ever published of X- 
ray measurements of this class of objects, the sample is 
small and admittedly far from complete or unbiased. On 
the contrary, it is clearly biased toward X-ray bright- 
est (and/or therefore potentially less obscured) objects. 
Moreover, an accurate determination of the ionization 
state of the absorber is beyond the capability of the cur- 
rent high-resolution X-ray detectors for objects of X-ray 
flux ~ 10~ 13 erg cm -1 s , and must await the next gen- 
eration of X-ray missions. Questions such as the nature 
of the gaseous nuclear environment of the early evolution- 
ary stages of AGN Ijdi Matteo et al. 2005|) . or the feed- 
back between the (recurrent? |Siemiginowska et al. 2 002 
2003) AGN activity and its environment, or surrounding 
clusters ( Fa bian et al. 20031 |Siemig inowsh a et al. 20 05a) 
could be properly addressed by high quality spectroscopic 
X-ray measurements of a sizable sample of GPS sources. 
Nonetheless, the results presented in this paper open a so 
far unexplored window on the environmental conditions 
surrounding the origin of the radio power in the Universe. 
An extension of this study to a complete and unbiased 
sample of GPS galaxies is under way. 
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